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RESEARCH MEMORANDUM 


THE EFFECTS OF FLEXIBILITY ON THE LONGITUDINAL AND LATERAL- 

DIRECTIONAL RESPONSE OF A LARGE AIRPLANE 

By Henry A. Cole, Jr,, Stuart C. Brown, 
and EJuclld C. Holleman 


SUMMARY 


Longitudinal and lateral-directional frequency responses determined 
from transient fli^t data excited "by control pulses are presented for 
a large flexible swept-wing airplane for a wide range of flight condi- 
tions. Transient flight data are reduced to transfer-function form and 
comparisons are made with predicted values for the rigid and flexible 
airplane. Flexible -airplane transfer-function coefficients detezmined 
by the pseudo static method show good agreement with those dete rmin ed 
from transient flight data. 

The effects of structural modes on the frequency response at various 
locations on the airplane are shown. These responses are compared with 
node lines of the calculated normal free -free modes of the airplane with 
particular reference to optimum location for automatic contiml pickups. 
Good correspondence is obtained between the measured frequency responses 
and the node lines calculated by dynamical analyses. 


INTRODUCTION 


In the design of an automatic control system for an airplane, the 
airplane response to control-surface motions should be known accurately. 
In the past the response has been fairly well established for relatively 
rigid airplanes by flight tests and theory (refs. 1 and 2). However, in 
recent years, the desire to increase the range and speed of large air- 
planes has led to sweptback wings of high aspect ratio, thin airfoils, 
and fuselages of high fineness ratio. All of these factors tend to 
increase the flexibility of the structure and the associated aeroelastic 
effects are becoming of greater importance in problems of static and 
dynamic stability and control. The dynamic effects are especially impor- 
tant in the design of automatic control systems because structural modes 
may introduce instabilities which would not axlse with a rigid airplane. 
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Hence, It Is Important for the automatic control designer to consider 
the effects of flexihlllty on his control systems . 

NOTATION 

Cj^ slope of lift curve 

D differential operator, ~ 

Ki rolling-velocity gain 

Kg angle -of -yaw gain 

Ky radius of gyration about principal lateral axis, mean aerodynamic 

chords 

Kg pitching-velocity gain 

M Mach ntimber 

Tg pitching-velocity time constant, sec 

n normal acceleration, positive downward, gravity units 

p rolling velocity, radlans/sec 

q dynamic pressure, Ib/sq ft 

r yawing velocity, radlans/sec 

t time 

6a total aileron deflection, radians 

&e elevator deflection, radians 

6j. rudder deflection, radians 

6 pitching velocity, radians/sec 

^ damping ratio, dimensionless 

characteristic numerator constant, dimensionless 
cj frequency, radlans/sec 

£ 0 n undamped natural frequency, radlans/sec 
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(jJa characteristic numerator constant, radians/sec 

T rolling-velocity time constant, sec 


DISCUSSION 


In the synthesis of an automatic control system, the aerodynamics 
are usually represented "by a transfer function in a "block diagram. As 
an example, the transfer functions of pitching velocity to elevator 
control for a rigid and a flexi"ble airplane are presented in figure 1. 

The complications introduced "by flexi"bility are apparent. Not only are 
the coefficients of the similar teims different from those of the rigid 
airplane, "but also a large number of second -order terms are added to the 
numerator and denominator by the structural modes. These terms cannot 
be taken li^tly because they have been directly responsible for a num- 
ber of instabilities in control systems in the past (refs. 3 and It). 

Also, analysis of autopilots in vhlch control servos introduce similar 
terms has shown that the second-order denominator terms limit the allow- 
able autopilot gain for stability. Fortunately, only the terms which 
are significant within the frequency range of the control servo need to 
be included. Some frequency-response data on the dynamic characteristics 
of the control system of a typical large airplane which may be useful in 
autopilot design are given in reference 5. 

In order to enlarge on the experience and analytical techniques 
necessary to evaluate these aeroelastic effects, the NACA is fll^t test- 
ing a Boeing airplane over its operating range of altitudes and 

Mach numbers. The entire structure was Instrumented with strain gages, 
accelerometers, and an optigraph as well as standard instruments so that 
the response of the complete airplane could be measTired. Some of the 
measured acceleration frequency responses to elevator motion are shown 
in figure 2 . The amplitude ratio is given on the upper plot and the phase 
angle on the lower one for a range of frequencies including the short- 
period mode, from 1 to 3 radians per second, and a number of structural 
modes, from 8 to 30 radians per second.^ The acceleration responses at 
the center of gravity, wing tip, tail, and a point on the wing near the 
Inboard nacelle are given. Several features of the response of the test 
airplane are apparent from these curves: (l) At frequencies below 

6 radians per second, the neighborhood of the short period mode, the 
response is similar for all locations on the airplane and (2) at frequen- 
cies above 6 radians per second, the region of the structural modes, 
the response differs widely with location on the aizqalsine and la char- 
acterized by multiple peaks in amplitude ratio and large shifts in phase 

^Note that the amplitude ratio and frequency are plotted to loga- 
rithmic scales so that the frequency-response curves can be Interpreted 
in terms of the characteristic forms of the transfer-function terms 
given in reference 6. 
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angles. The peaks correspond to the second -order denominator tenna and 
the valleys to the second -order numerator terms in the flexible-airplane 
transfer fimctlons. The part of the response below the structural modes, 
which is of most concern to us from a maneuvering and control standpoint, 
will be discussed first and the region of the structural modes later on. 


Low-Frequency Response 


The question is: How well can this low-frequency response be pre- 
dicted? One way is to use the same form of transfer functions as for 
the rigid airplane, but to modify the coefficients for flexibility effects. 

In this method, which is called pseudostatic, the individual stability 
derivatives which are used to calculate the transfer-function coefficients 
are determined for the condition of static equilibritim between the struc- 
tural spring forces and the aerodynamic loads. In addition to the usual 
derivatives of rigid airplanes, the lifts and moments resulting from 
deflections due to inertial loadings (normal acceleration and pitching 
acceleration) are Included. These derivatives resulting from inertial 
loads are very important and vary widely depending on the mass distribu- 
tion of the airplane. The pseudostatic method and its application to 
this airplane is described in detail in reference 7 - Extensive work on 
aeroelastic effects on this airplane has also been reported in reference 8. 
Pitching- velocity responses predicted fon pseudostatic, and rigid airplanes .. 

are compared with the response measured in fli^t in figure 3* The agree- 
ment of experiment with the pseudostatic method is very good up to about 
5 radians per second after which the structural modes dominate the s 

response. The predicted rigid-airplane response, falls co.nslderably 
below the experimental response. It should be realized that accurate 
predictions of these responses depend a good deal on accurate estimates 
of the spring and mass characteristics of the airplane. In the analysis 
used here, the spring characteristics of the airplane were determined 
from static deflection measurements reported in reference 9, and the mass 
characteristics by ground oscillation tests which are as yet unpublished. 


Longitudinal Transfer -Function Coefficients 


A more comprehensive picture of the effects of flexibility over the 
fli^t range is obtained by examining the transfer-function coefficients. 

In figure V, the denominator transfer-function coefficients of figure 1 

are plotted versus the ae.roelastlc parameter, dynamic pressure multiplied 

by the ratio of the slopes of lift curve at the test Mach number and at 

zero Mach nirnber. Also curves for 15,000 and 35>000 feet are shown so 

that the flight range of the airplane is pretty well covered. The solid « 

line represents the predicted pseudostatic values and the dashed curve 
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■the predicted rigid ones. Experimen'tal points sTaown on "these curves were 
determlneil from the free oscillations foUo'vring pulse-elevator Inputs. 

The predicted values Indicate a large reduction in natural frequency, 
tiijx, due to flexibility and this trend is veil supported by the experi- 
mental points. This reduction in natural frequency is principally a 
result of the loss in the spring term, due to fuselage and stabilizer 

bending. The loss in frequency would be even greater if it were not for 
compensating pitching moments contributed by ■wing bending due to normal 
acceleration, a result primarily due to placing the nacelle masses near 
the wing tip. 

As for the damping ratio, the damping and spring terms are reduced 
by flexibility but the compen sating effects of inertial loads are such 
that there is little difference between the predicted rigid and pseudo- 
static values. The level of the experimental points is in good agreement 
with the predicted ones, but there is considerable scatter, probably 
because measurements of "this coefficient are ill-conditioned when the 
transient motion damps out in a couple of cycles or less. 

The numerator coefficients are shown in figure 5- S®iHj Kg, is 

mainly dependent on "the ratio of the forcing term Cmg, to the spring 
term for the rigid airplane . In the flexible airplane there are 

losses in bo'th "the forcing term and the spring "teim due to structural 
bending as well as compensating effects of inertia as mehtioned in con- 
nection with the natural frequency. Hence, althou^ there are large 
changes in the indi"vidual derivatives due 'to flexibility, it Just happens 
in this case "bhat the over-all effect on the gain is small. The experi- 
mental points were obtained from slow push-pull maneuvers and they agree 
well with the predicted ones. 

The time constsint, Tg, is primarily the ratio of "the mass to the 
damping force in translation and represents the phase lead of pitching 
velocity to angle of attack. This term is increased by flexibility 
because of the reduction in the damping force resulting from wing bend- 
ing. The experimental points were obtained by the pulse technique and 
agreement with the pseudos"tatic values is good at 15,000 feet. At 
35,000 feet the agreement is not very good which is probably due to Mach 
nimiber effects ■which are more impor"tant than aeroelastic effects at this 
altitude . 


Lateral Response 


Similar analyses and flight-test measurements are being made for the 
lateral -directional response to aileron and rudder controls. Typical 
frequency responses (yaw damper off) obtained by the pulse technique are 
shown in figures 6 and 7 - Also shovm on "these curves for comparison are 
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the predicted pseudostatlc and rigid -ali^p lane responses. Figure 6 pre- 
sents the frequency response of rolling velocity to aileren control and 
figure 7 "the yawing-velocity response to rudder control. The first peak 
in either figure represents the Dutch ivoU mode, and the peaks above 
5 radians per second the structural modes. The aileron control excites 
the wing first-antisyimiietrical mode primarily, and the rudder control 
excites the hody side-hending mode. As in the longitudinal case, the 
response at low frequencies is similar to that of the rigid airplane for 
the structural modes are well separated from the Dutch roll mode. Hence, 
the pseudostatic approach should also apply to the low-frequency range 
of the lateral-directional response. 


Latereil -Directional Transfer-Function Coefficients 


Some preliminary values of transfer-function coefficients have been 
obtained by matching transient time histories from control pulses on sin 
electronic analog. The transfer-function; coefficients used for this 
work are shown in figure 8, and are the simplified transfer functions 
used for rigid airplanes. The denominator coefficients are shown in 
figure 9. The lateral-directional fli^t testing has not been completed 
so data are only shown for an altitude of 25,000 feet. First, consider 
the Dutch roll mode which is represented by 'ttie second-order denominator 
term in the transfer functions of figure 8. The frequency, tiJn, is lowered 
primarily by a loss in the spring term, Cnp, due to fuselage and fin 
bending. The damping ratio, is also lowered primarily by fuselage and 
fin bending which results in a loss in damping in yaw. The time constant, 
T, depends primarily on the ratio of the moment of inertia in roll to the 
damping in roll. While damping in roll is reduced by wing bending, the 
effective inertia in roll is also reduced by the rolling moments Induced 
by wing bending under rolling acceleration loadsj; hence^ the over-all 
effect is small. The experimental points show good agreement with the 
predicted pseudostatic values. Note that there is much less scatter in 
the damping ratio than for the longitudinal case because we have many 
cycles in the transient time histories. 

The numerator terms are plotted in figure 10. The rolling-velocity 
gain, Kx, is decreased considerably by wing bending. Predictions indi- 
cate that this gain goes through zero at a dynamic pressure of about 
720 pounds per square foot. This point, of course, is known as aileron 
reversal and will be approached more closely in tests at lower altitudes. 
The aileron reversal is most critical for outboard ailerons, which were 
used on this airplane (ref, lO) , and for this reason there has been a 
trend to inboard ailerons and spoilers for lateral control of high-speed 
airplanes. -- - - - ._l. . 

It may he shown that, when the flight-path angle is zero and the 
angle of attack is low, the terms o)a same as Uq and 
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respectively, to a first-order approximation (ref. 2 ) . By comparison of 
values of Wa and with those of and ^ in figure 9, it may be 
seen that they are nearly the same when is large. Hence, the second- 
order numerator term nearly cancels the Dutch roll term, and the rolling- 
velocity transfer function is sometimes simplified to Ki/(l + tD) • 

The numerator constant. Kg, in the yawing-velocity transfer function 
Is primarily the forcing term, , over the spring term. Cup. The tail 

contribution of each is reduced about the same due to fuselage and fin 
bending, but Cno also has a fuselage contribution, so the over-all effect 
of flexibility ii to increase K^. From these preliminary data, it appears 
that the pseudostatic method is also adequate for the prediction of low- 
frequency response for the lateral-directional case. 


Experimental Technique 


Although it is not the Intention of this paper to give detailed 
information on the reduction of data, the reasons for the different 
methods used in determining experimental transfer- function coefficients 
should be of interest. Generally speaking, a large amount of scatter 
results in some of the coefficients when they are all determined from 
one transient record because the conditioning of the coefficients varies 
with the frequency content of the input. Greater accuracy is obtained 
if the ’coefficients are determined from several inputs which cover the 
frequency range of interest. Hence, gain constants should be determined 
from slow maneuvers. Denominator terms should be determined from free 
oscillatlona following inputs which excite primarily one mode at a time. 
When accurate values of gain and denominator coefficients are known, then 
the other nimierator terras can be determined from the forced oscillation 
part of pulse control inputs. This work can either be done on an 
electronic analog or by curve-fitting the frequency response. 


Hi^-Frequency Response 


The hi. ^-frequency part will now be examined in more detail, by 
referring again to figure 2. An acceleremeter pickup in an automatic 
control system, if located at the wing tip, center of gravity, or tail, 
would feed voltages into the system in accordance with the hi^ peaks and 
large phase shifts in the frequency responses. As mentioned before, this 
would be apt to compromiee the stability of the system if the control 
servo could respond to these frequencies. However, we note that if the 
accelerometer were placed on the wing near the inboard nacelle where the 
structural -mode peaks are absent, the stability of the system would not 
be compromised. Optimum locations for pickups such as this can be deter- 
mined from fli^t tests, but the question arises: How well can they be 

determined prior to fli^t? 
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In order to check this, the normal free -free modes of the complete 
airplane were calculated by methods similar to those in reference 11. 

The node lines for the first three modes are presented in figure 11. The 
first mode at a frequency of 8 radians per second is primarily wing first- 
bending mode and is called "wing flapping" by the pilots. In figure 2, 
it is noted that the response only peaks up near- the wing tips at this 
frequency and that the response is small for locations on the fuselage near 
the node lines. The next predicted node at 23 radians per second is pri- 
marily wing first torsion. Then at a natural frequency of 27 radians per 
second is the body first-bending mode. It is noted that a large double 
peak occurs on the measured -frequency responses in the range from 20 to 
30 radians per second which is undoubtedly due to the occurrence of these 
two modes at nearly the same frequency. Also, the peak is absent at the 
point where the node lines cross on the wing near the Inboard nacelle. 
Hence, there appears to be good correspondence between the flight node 
lines and those predicted by dynamic analysis. 

There are several other peaks in the wing-tip acceleration response 
besides those mentioned. Some of these are believed to be due to flexi- 
bilities in the nacelle mountings which were neglected in the analysis. 
Also, an additional mode appears near I 6 radians per second when the 
short-period mode is included in the analysis (see ref. 7 ) • 

This means that the effects of flexibility can be taken entirely 
into account by pseudostatic methods if pickups are located on these node 
lines. If this is not possible, then the structural mode frequencies 
would either have to be filtered out or the control system tailored to be 
stable with their presence. It certainly would not be practical to locate 
all pickups and instruments at this point near the nacelle so the next • 
beat thing would be to decide which mode would affect the control system 
the most and to locate the pickup on the node line for that mode. For 
example, the high peak at 20 to 3 O radians per second could also be 
avoided where the node lines cross the fuselage near the tall. The lines 
shown here are for accelerometer pickups, but lines of zero angular- 
velocity response can also be determined from the calculated modes for 
use in locating turnmeter pickups. 

COHCmDING REMAEKS 


In general, the pseudostatic method adequately predicts the longi- 
tudinal and latersl-dlrectional frequency response in the region below 
the structural -mode frequencies. Althougji. in some cases there were 
only small apparent differences between rigid and pseudostatic transfer- 
function coefficients, the individual aeroelastlc effects on the wing 
and tall are large and depend a good deal on the airplane mass distri- 
bution. The agreement of experiment with pseudostatic predictions 
Indicates that these large compensating effects are properly evaluated 
by the pseudo static method. 
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The correspondence of "bhe measured modes in fli^t; vith "those calcu- 
lated by dynamical analyses indicates that optinnm locations for control 
system pickups can be determined prior to fli^t. Vfhen pickups are not 
located on these node lines, the effects of the additional second -order 
terms in the flexible-aii^plane transfer functions should be taken into 
accoTont . 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Apr. l4, 1955 
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Figure 3 .- Pitching-velocity frequency response to elevator at c.g.; 
hp = 25,000 feet, M = 0.7, c.g. = 20-percent M.A.C.j, W = 108,000 
pounds. . .. • 
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Figure 4.- Denominator coefficients of pitching-velocity transfer func 
tionj c.g. = 21-percent M.A.C., W =115,000 pounds, Ky^ = 2.4. 
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Figure 5 
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EXPERIMENT 



Figure 6.- Rolling -velocity frequency response to aileron at c.g.; 
hp = 35,000 feet, M = 0.6, W = 117,000 pounds. 
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Figure 7>- Yawing-velocity frequency response to rudder at c.g.; 
hp = 24,000 feet, M = 0.64, W = 104,000 pounds. 
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Figure 8.- Lateral-directional, trsinsfer functions. 
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Figure 9*~ Denominator coefficients of lateral-directional transfer 
functions; lip - 25,000 feet, W = 115,000 poimds. 
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Figure 10.- Numerator coefficients of lateral-directional transfer func- 
tions; hp = 25,000 feet, W = 115,000 pounds. 
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Figure 11.- Calculated longitudinal node lines. 
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